﻿// Histogram.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。
//

#include "pch.h"
#include <iostream>
#include "opencv2/opencv.hpp"
using namespace cv;
class Histogram1D {//一个柱状图
private:
	int histSize[1];//高度
	float hranges[2];//高度范围
	const float* ranges[1];//范围
	int channels[1];//通道
public:
	Histogram1D() {//初始化
		histSize[0] = 256;//灰度值大小为0~255总共256个
		hranges[0] = 0.0f;
		hranges[1] = 255.0f;
		ranges[0] = hranges;
		channels[0] = 0;
	}
	MatND  getHistogram(const Mat &image) {
		MatND hist;
		calcHist(&image, 1, channels, Mat(), hist, 1, histSize, ranges);
		return hist;
	}
	Mat getHistogramImage(const Mat &image) {
		MatND hist = getHistogram(image);
		double maxValue = 0.0f;
		double minValue = 0.0f;
		minMaxLoc(hist, &minValue, &maxValue);
		Mat histImg(histSize[0], histSize[0], CV_8U, Scalar(255));
		int hpt = static_cast<int>(0.9*histSize[0]);
		for (int h = 0; h < histSize[0]; h++) {
			float binValue = hist.at<float>(h);
			int intensity = static_cast<int>(binValue*hpt / maxValue);
			line(histImg, Point(h, histSize[0]), Point(h, histSize[0] - intensity), Scalar::all(0));
		}
		return histImg;
	}
};
int main()
{
	char fn[] = "D:\\opencv\\sources\\samples\\data\\pic2.png";
	Mat img = imread(fn,0);
	imshow("pic2", img);
	Histogram1D hist;
	Mat histImage;
	histImage = hist.getHistogramImage(img);
	imshow("histo", histImage);
	Mat bw;
	threshold(img, bw, 0,255,THRESH_OTSU);
	imshow("bw", bw);
	waitKey();
}

// 运行程序: Ctrl + F5 或调试 >“开始执行(不调试)”菜单
// 调试程序: F5 或调试 >“开始调试”菜单

// 入门提示: 
//   1. 使用解决方案资源管理器窗口添加/管理文件
//   2. 使用团队资源管理器窗口连接到源代码管理
//   3. 使用输出窗口查看生成输出和其他消息
//   4. 使用错误列表窗口查看错误
//   5. 转到“项目”>“添加新项”以创建新的代码文件，或转到“项目”>“添加现有项”以将现有代码文件添加到项目
//   6. 将来，若要再次打开此项目，请转到“文件”>“打开”>“项目”并选择 .sln 文件
